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ABSTRACT 

We report the discovery of a young massive stellar cluster and infrared nebula in the direction 
of the CS molecular cloud associated to the IRAS point source 16132-5039. The analysis of the 
mid-infrared images from the more accurate MSX catalog, reveled that there are two independent 
components associated with the IRAS source. The integral of the spectral energy distribution for 
these components, between 8.28 fim and 100 /xm, gave lower limits for the bolometric luminosity 
of the embedded objects of 8.7 x lO^L© and 9 x 10"^ Lq, which corresponds to ZAMS 08 and 
BO. 5 stars, respectively. The number of Lyman continuum photons expected from the stars that 
lie along the reddening hue for early-type stars is about 1.7 x 10^^ s~^, enough to produce 
the detected flux densities at 5 GHz. The NIR spectrum of the nebula increases with frequency, 
implying that free- free emission cannot be the main source of the extended luminosity, from which 
we conclude that the observed emission must be mainly dust scattered light. A comparison of 
the cluster described in this paper with the young stellar cluster associated with the IRAS source 
16177-5018, which is located at the same distance and direction, shows that the mean visual 
absorption of the newly discovered cluster is about 10 magnitudes smaller and it contains less 
massive stars, suggesting that it was formed from a less massive molecular cloud. 



Subject headings: stars 
ISM: dust, extintion 

Introduction 
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Massive stars are born within dense molecular 
clouds forming clusters and associations; during 
their formation and early evolution they are of- 
ten visible only at infrared and radio wavelengths. 
With the advent of the large bidimensional near- 
infrared array detectors, the morphological and 
photometric studies of extremely young galactic 
stellar clusters have benefited. At near-infrared 
wavelengths (1 to 2.5 /im) it is possible to probe 
deep into the dense dust clouds where star forma- 
tion is taking place. 

The strong ultraviolet (UV) radiation emitted 
by the massive stars dissociates and ionizes the 



^ Based on observations made at Laboratorio Nacional 
de Astrofisica/MCT, Brazil 



gas, forming compact HII regions seen at radio 
wavelengths. A large fraction of the radiation 
could also heat the dust, which eventually radiates 
in the far-infrared (FIR); for this reason, compact 
HII regions are among the brightest and most lu- 
minous objects in the Galaxy at 100/im. 

Because massive stars evolve very fast, they are 
extremely rare and difficult to find, except when 
related to the emission of the surrounding molecu- 
lar cloud. In that sense, CS and NH3 lines at radio 
frequencies, characteristic of high density gas, are 
good tracers of massive star forming regions. 

In this paper we present observations in the 
near infrared, of a young stellar cluster of mas- 
sive stars in the direction of the IRAS point 
source 116132-5039. This work is a part of a 
survey aimed to the identification of stellar pop- 
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ulations in the direction of IRAS sources that 
have colors characteristics of ultracompact HII 
regions (Wood & Churchweh 1989) and strong 
CS (2-1) Hne emission (Bronfman ct al. 1996). 
The studied region is located in the direction 
of another massive star forming region asso- 
ciated to the IRAS point source 16177-5018 
(Roman-lopes et al. 2003); they are part of the 
RCW 106 complex, located in the southern Galac- 
tic plane at a distance of 3.7 kpc (Caswell & Haynes 1987) 
The near-infrared cluster presented here had re- 
cently been located by Dutra et al. (2003) by 
visual inspection of the 2MASS images; however 
their work contains only the possible cluster loca- 
tion, and as Persi, Tapia & Roth (2000) already 
showed for NGC6334IV, a high concentration of 
K band sources can be due to localized low ex- 
tinction and be mistaken with a stellar cluster, 
leading to false identifications. 

The cluster is associated with the radio source 
332.541-0.111, which presents continuum radio 
emission at 5 GHz as well as hydrogen recombina- 
tion lines (Caswell & Haynes 1987). The observa- 
tions and data reduction are described in section 
2, the results are presented in section 3 and our 
main conclusions are sumarized in section 4. 

2. Observations and data reduction 

The imaging observations were performed in 

June 2001 and May 2003 with the Near Infrared 
Camera (CamlV) of Laboratorio Nacional de As- 
trofisica (LNA), Brazil, equipped with a Hawaii 
1024x1024 pixel HgCdTe array detector mounted 
on the 0.6 m Boiler & Chivens telescope. The ob- 
servations consisted of 8'x 8' frames in the direc- 
tion of the IRAS source 16132-5039; the plate scale 
was 0.47 arcsec/pixel and the mean values of the 
PSF full width at half maximum (FWHM) were 
1.4, 1.7 and 2.1 arcsec at the J, H and nhK im- 
ages. The total integration time was 2220 s for 
J, 1440 s for H and 6400 s for nhK filters, result- 
ing in a sensitivity at 3(t of 18.2, 17.4 and 14.2 
magnitudes and completeness limits of 17.8, 16.5 
and 13.8 magnitudes respectively. Details about 
the calibration and reduction procedures can be 
found in Roman-Lopes et al. (2003). 

Photometry from 2MASS All Sky Point Source 



Catalogue ^ in the J, H and Ks filters became 
available recently, with completeness limits of 
15.8, 15.1 and 14.3 magnitudes, in the three pass- 
bands (Egan et al. 2001). Since the 2MASS Kg 
band photometry, centered at 2.17 /um and with 
a bandpass of 0.32 /zm has a completeness limit 
greater than our libK photometry, we completed 
our observations with the 2MASS catalogue as- 
trometry and magnitudes in this filter. The com- 
parison of our photometry with the 2MASS survey 
data, in an area of about 40 square arcmin, corre- 
sponding to a total of 469, 610 and 440 common 
sources in the J, H and K bands, respectively, is 
shown in figure 1. We see a good linear relation 
between the two systems, with a slope of 1 and a 
dispersion that increases with the magnitude. 

3. Results &: Discussion 

The combined false-colour infrared image, (J 
blue, H green and nbif red) of the whole field is 
displayed in Figure 2, together with amplified in- 
dividual images of the nebular region at all bands. 
All of them, but especially the H and nbif images, 
shows the presence of a small spheroidal nebula 
with a bright star at its center. 

3.1. The IRAS source 

In Figure 3 we present a contour map con- 
structed from the LNA rHoK image, with a beam 
size of 2 X 2 pixels which shows the region around 
IRAS16132-5039. The contours start at 2.2x10"'' 
Jy/beam and are spaced by this same value. The 
IRAS coordinate has an intrinsic error delimited 
by the ellipse plotted in the figure. A more accu- 
rate position for the IR source was obtained from 
the Midcourse Space Experiment - MSX point 
source catalog {psc) ^. The MSX surveyed the 
entire Galactic plane within | 6 | < 5° in four mid- 
infrared spectral bands centered at 8.28, 12.13, 
14.65 and 21.34 /zm, with image resolution of 
19 arcsec and a global absolute astrometric ac- 
curacy of about 1.9 arcsec (Price et al. 2001). 
We found one MSX source, with coordinates 
a(J2000) = 16^17™02.47^ <5(J2000) = -50<i47'"03.5^ 
within the IRAS uncertainty ellipse; its coordi- 
nates coincides with the star we labeled IRSl. 



^http:/ /www. ipac.caltcch.edu/ 

^http;/ /www. ipac.caltech.edu/ipac/msx/msx. html 
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However, looking at the MSX images, we 
found another closeby source, although outside 
the IRAS uncertainty ellipse, with coordinates 
a(J2000) 16'^16™55.94% (5(J2000) = -50'^47'"07.8^ 
In Figure 4 we present our H band image over- 
laid with the 8.28/xm band MSX contour diagram, 
with the contours spaced by 8 x 10^^ W m~^ 
sr-\ starting at 2.8 x IQ-^ W m'^ sr"!. The 
same source was found in the contour plots from 
the other MSX bands. We designated the stronger 
source as A and the other as B. While the IRSl 
object and the infrared nebulae are related to 
source A , source B is also associated with a small 
nebular region that shows a concentration of NIR 
sources, as can also be seen in Figure 4. 

Since only the 8.28 fim flux density was given 
in the MSX psc for both sources, we integrated 
the flux density of each individual source for the 
four mid-infrared bands; the results arc presented 
in table 1, which also shows the values from MSX 
and IRAS catalogs. Our integrated flux density 
for source B in the 8.28 /xm coincides within 5% 
with the value given in the MSX catalogue, but 
it is 30% larger for the A source. This discrep- 
ancy can be understood if we consider that the 
automatic MSX algorithm subtracted part of the 
source B flux density as background contribution. 
The relative contribution of source B dcxTcases 
with increasing wavelength, explaining also why it 
was not resolved by the MSX algorithm. It should 
be noticed that the reported IRAS flux density at 
12 /xm coincides with the sum of our derived flux 
densities at 12.13 /zm within 10%, while the MSX 
value was about 50% lower. 

In Figure 5 we plotted the mid to far-infrared 
spectral energy distribution of the sources A and 
B, without any correction for absorption. We as- 
sumed that the IRAS flux density in the far in- 
frared is divided between the two sources in the 
same way as in the mid-infrared (14.65 and 21.33 
Ilia): about 90% originating from source A and 
10% from source B. We then integrated the ob- 
served flux densities in the mid-far infrared, as- 
suming a distance of 3.7 kpc, obtaining a lumi- 
nosity La = 8.7 X IO^'Lq and Lb = 9 x IO^Lq 
for sources A and B, respectively. Assuming that 
the IR flux density represents a lower limit for 
bolometric luminosity of the embedded stars, we 
derived ZAMS spectral types 08 and BO. 5 for 
the sources A and B, respectively (Hanson et al. 



1997). 

3.2. Cluster population 

In order to examine the nature of the stellar 

population in the direction of the IRAS source, we 
analyzed the stars in two delimited regions: one 
that we labeled "cluster" , which contains the neb- 
ula and another that we labeled " control" , which 
has a stellar population that we believe is domi- 
nated by "field" stars, as illustrated in Figure 6. 
We represented the position of all objects detected 
in the H band by crosses and we can see that the 
small region labeled "cluster" shows a concentra- 
tion of sources. In figure 7 we present comparative 
(J — H) versus {H — K) diagrams for the stars de- 
tected in the J, H and K images, together with 
the position of the main sequence, giant branch 
and reddening vectors for early and late type stars 
(Koornneef 1983). 

We see that the stellar population of the "clus- 
ter" region is quite different from that of the " con- 
trol" region, with many sources lying on the right 
side of the reddening vector for early type stars, 
showing excess emission at 2.2 /xm. In the "con- 
trol" region the majority of the sources have col- 
ors of reddened photospheres, with many objects 
located along the reddening vector for late type 
stars. It is well established that very young pre- 
main sequence objects present large infrared ex- 
cess due to the presence of warm circumstellar 
dust (Lada & Adams 1992). Our results suggest 
that the stellar population in the direction of the 
IRAS source is very young, as can be inferred from 
their position in the (J — H) versus [H — K) dia- 
gram. 

We separated the cluster sources from the field 
stars, by selecting all sources that lie to the right or 
on the reddening vector for early type stars in the 
cluster's color-color diagram. Table 2 shows the 
coordinates and photometry of all selected sources 
(J, H and K magnitudes from both LNA and 
2MASS surveys). 

Further information about the nature of the se- 
lected objects in Table 2 can be extracted from 
J versus (J — H) color-magnitude diagram shown 
in Figure 8. We used this diagram instead the K 
versus [H — K) color-magnitude diagram to min- 
imize the cfcct of the "excess" of emission in the 
NIR on the derived stellar spectral types. The lo- 
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cus of the main-scqucncc for class V stars at 3.7 
kpc (Caswell & Haynes 1987) is also plotted, with 
the position of each spectral type earlier than AO 
V indicated. The intrinsic colors were taken from 
Koornneef (1983) while the absolute J magnitudes 
were calculated from the absolute visual luminos- 
ity for ZAMS taken from Hanson et al. (1997). 
The reddening vector for a ZAMS BO V star, taken 
from Rieke & Lebofsky (1985), is shown by the 
dashed line with the positions of visual extinc- 
tions Ay = 10 and 20 magnitudes indicated by 
filled circles. We also indicated the sources with 
and without "excess" in the color-color diagram 
by open and filled triangles respectivelly. 

For the assumed distance, we estimated the 
spectral type of the stars that do not present ex- 
cess emission in the near infrared, by following the 
de-reddening vector in the color-magnitude dia- 
gram, for the others we only gave a rough classi- 
fication; the results are shown in the last column 
of Table 2. The main source of error in the de- 
rived spectral types arises from uncertainties in 
the assumed cluster distance, which was derived 
from the velocities of the radio hydrogen recombi- 
nation lines. Since the closest distance was used, 
the main uncertainties come from the errors in the 
galactic rotation curve model and related param- 
eters. From the works of Blum et al. (1999, 2000, 
2001) and Figueredo et al. (2002) comparing kine- 
matic with spectroscopic distances we find that 
they do not differ in more than 1 kpc, in which 
case the change in the luminosity class would be 
of two sub-spectral types for early O stars and one 
for early B stars. 

We must notice that source IRSl, which is as- 
sociated with MSX source A, has an estimated 
spectral type of at least 05, reddened by about 
Ay ~ 14 magnitudes. Besides, there are five ob- 
jects (IRS3, IRSll, IRS18, IRS21 and IRS33) as- 
sociated with the mid-infrared source B; IRS3 is 
probably an 08 ZAMS star reddened by about 
Av= 7 magnitudes, while the others have esti- 
mated spectral types of early-B stars. The corre- 
sponding bolometric luminosities are in agreement 
with the lower limits derived from the integrated 
mid-far infrared flux densities, corresponding to 
08 and BO. 5, for sources A and B respectively, as 
seen in section 3.1. 

A lower limit to the number of Lyman-continuum 
photons produced in the star forming region, can 



be calculated taking into account only the stars 
that do not show "excess" in the color-magnitude 
diagram (IRSl, 6, 8, 9, 10, 13, 14, 16, 20, 23, 26, 
27 and 35 in Table 2). It was computed from the 
relation given by Hanson et al. (1997), resulting 
in 1.7 X 10^® photons s~^. Is interesting to note 
that IRSl is responsible for more than 90% of 
the Lyman continuum photons, being the main 
ionization source in the whole region. 

It is also possible to obtain the number of ion- 
izing Lyman continuum photons N^y from the ra- 
dio continuum flux density given by Caswell & 
Haynes (1987), using the expression derived by 
Rubin (1968): 

_ 5.59 X 104Sg(^) j^2y-0.45^0.1 

l + U[<He+l(H+ + He+)] ^ > 

where v is in units of 5 GHz and /j is the fraction of 
helium recombination photons that are energetic 
enough to ionize hydrogen. Using the values of 
Te = 4500 K, S{v) = 3.3 Jy and D= 3.7 kpc taken 
from Caswell & Haynes (1987), we find that « 
and NLy 6 x 10^^ photons s^^, compatible with 
the lower limit derived from the observed stars. 

3.3. The Infrared Nebula 

We can see from the; detailed J, H and vibK 
images in figure 2 that the nebula presents a 
spheroidal shape, approximately symetric around 
the IRSl source. In figure 9 we shown contour 
maps of the nebular region at J, H and nbi^ bands 
obtained from our infrared images. The contours 
were calibrated in flux with the values starting 
at 0.44 (J), 1.6 {H) and 2.2x10-" Jy/beam {K), 
with intervals of 0.37, 1.8 and 1.1x10"" Jy/beam 
respectively. We estimated the total flux density 
by measuring the area between contours and mul- 
tiplying by the value of the corresponding flux 
density per imit area, obtaining S{J) = 0.023 Jy, 
S{H) = 0.07 Jy and S{K) = 0.11 Jy. We then cor- 
rected these results for extinction, using the mean 
value < Ay > = 15.1 magnitudes, derived from 
the stars in the direction of the nebula that do not 
present infrared excess (IRSl, IRS6 and IRSIO), 
and the standard extinction law from Rieke & 
Lebofski (1985). We obtained 5(J) = 1.16 Jy, 
S{H) = 0.8 Jy and SiK) = 0.52 Jy. 

In the previous section we showed that the num- 
ber of ionizing photons available from the detected 
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stars is enough to explain the radio continuum flux 
density measured by Caswell & Haynes (1987). 
We will investigate now the contribution of free- 
free emission to the observed nebular IR flux den- 
sity. Assuming constant density and temperature 
across the cloud and local thermodynamic equilib- 
rium, the flux density S{i') due to free-free emis- 
sion can be written as: 



Siu) = T,B,{T)n 



(2) 



where is the optical depth at frequency v, 
B^iT) is the Planck function 



2hv^ 



1 



c2 exp(/ij//A:T) - 1 
and n is the solid angle of the source given by: 
O = 7r(i/2£))^ 



(3) 



(4) 



where L is the diameter of the ionized cloud and 
D the distance to the observer. The optical depth 
at a given frequency v is: 



Tv = cxgL 



(5) 



where ag is the free-free absorption coefficient 
(cgs) taken from Rybick (1979): 



as 



3.7 X 10^ [1 - exp{-hi^/kT)] UeUi gs{i',T) 

1/3 z-2 ri/2 

(6) 

where gs{v,T) is the Gaunt factor obtained from 
Karzas & Latter (1961). 

For two frequencies vi and V2 the ratio of the 
corresponding flux densities S{yi) and S'(f2) may 
be calculated from: 

^(^^)-e.MM-2-.i)AT]^^ (7) 



S{V2) 



For a flux density in the 5 GHz continuum of 
3.3 Jy and an electron temperature T,, = 4500 K, 
as given by Caswell & Haynes (1987) the expected 
flux densities at the J,H and K bands are 0.05, 
0.10 and 0.16 Jy, respectively. We verify that the 
measured values are much larger than what was 
expected from the free-free emission, derived from 
the radio data. In fact, only an absorption as low 
as 4 magnitudes would explain the inferred spec- 
trum as free-free emission. Since in the direction 
of the A source, this absorption is incompatible 



with even the less absorbed star (IRSl), we be- 
lieve that this is not the main source of extended 
IR emission. 

Accepting the mean value of < Ay > = 15.1 
magnitudes as the mean value for the absorption 
in the direction of A source, we found that the 
corrected flux density increases with frequency, 
suggesting that the observed extended radiation 
is scattered light from the nearby stars. We ad- 
justed then a black body to the NIR fluxes, ob- 
taining a good fit for T« 16000K, characteristic 
of middle-B stars (Hanson et al. 1997). Lumsden 
& Puxley (1996), analyzing the ultracompact HII 
region G45. 12-1-0. 13, also obtained an extinction 
corrected flux density that increases with decreas- 
ing wavelengths and interpreted it as due to stellar 
light, scattered by dust through the HII region. 

3.4. Conclusions 

Near-IR imaging in the direction of the CS 
molecular cloud associated with the IRAS source 
16132-5039, revealed an embedded young massive 
stellar cluster. We detected 35 member candidates 
up to our completeness limit, concentrated in an 
area of about 2 square parsec. All images, but 
especially the H and nhK bands, show the pres- 
ence of a small spheroidal nebula with a bright 
star (IRSl) at its center. 

The stars associated with the IRAS point 
source were identified using more accurate po- 
sitions from the MSX catalogue. The analysis of 
the mid-infrared images reveled that there are two 
sources associated with IRAS 16132-5039. The 
strongest coincides with the position of at least 
a dozen of OB stars, while the weaker source is 
associated to less massive objects, with spectral 
types characteristic of middle-B ZAMS stars. The 
integral of the spectral energy distribution of the 
MSX-IRAS sources, between 8.28 /xm and 100 
fim, gives lower limits to the bolometric luminos- 
ity of the embeded objects of L = 8.7 x IO^'-Lq 
and L = 9 X 10^ Lq, which corresponds to ZAMS 
08 and BO. 5 stars, respectively (Hanson et al. 
1997). The results are compatible with the spec- 
tral types of the objects detected in the NIR, since 
it is possible that only part of the energy emitted 
by the stars is reprocessed by the dusty envelope. 
In that sense, they can be taken as lower limits to 
the bolometric luminosity of the embeded stars. 
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Assuming that the radio emission measured by 
Caswell & Haynes (1987) originates in this re- 
gion, at a distance of 3.7 kpc, we estimated the 
number of ionizing Lyman continuum photons as 
N^y ~ 6 X 10"'* photons s~^. On the other 
hand, the number of Lyman continuum photons 
expected from the stars that lie along the redden- 
ing line for early-type stars is about 1.7 x 10^^ 
s~^, enough to produce the detected flux densities 
at 5 GHz. The IRSl source is enough to account 
for more than 90% of the total number of Lyman 
continuum photons necessary to ionize the gas. 

Analysis of the integrated flux densities of the 
NIR nebula at the J, H and nbi^ bands revealed 
that they increase with frequency, implying that 
free- free emission cannot be the main source of the 
extended luminosity, unless we assume only four 
magnitudes of visual extinction. Since this value 
is incompatible with the extinction derived from 
the stars that do not shown excess of emission at 
2.2 yUTO, we conclude that the observed emission 
must be mainly dust scattered light. 

A comparison of the cluster described in this 
paper with the young stellar cluster associated 
with the IRAS source 16177-5018 (Roman-Lopes 
et al. 2003), which is located at the same dis- 
tance and direction, shows that the the former 
contains less massive stars. Since its mean visual 
absorption is also about 10 magnitudes smaller 
than that of IRAS 16177-5018, it is possible that 
it was formed from a less massive molecular cloud. 
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Table 1 

Integrated fluxes from MSX A, C, D and E images for the two mid-infrared sources (see 
text). Also are shown the data from MSX and IRAS point source catalogue {psc). 





Integration "A" 


Integration " B" 


MSX psc "A" 


MSX psc "B" 


IRAS psc 


8.28 


15.8 


4.3 


11.3 


4.1 




12 










44.8 


12.13 


32.8 


3.6 


23.0 






14.65 


23.9 


2.0 


21.3 






21.33 


117 


10.8 


123.8 






25 










266 


60 










2311 


100 










4618 
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Table 2 

List of the selected near-infrared sources 



IRS 


a(J2000) 


(5(J2000) 


JcamlV 




HcamlV 




KcamlV 




Spec Type 


1 


16:17:02.20 


-50:47:03.1 


12.26(4) 


11.97(4) 


10.88(4) 


10.65(4) 


10.07(4) 


9.73(4) 


05 


2 


16:17:09.23 


-50:47:14.7 


13.36(3) 


13.26(4) 


11.66(2) 


11.47(3) 


10.13(3) 


10.15(3) 


mid-O 


3 




— 50:47:23 


11.77(2) 




11.01(2) 




10.30(4) 






4 


16:17:00.62 


-50:47:49.6 


12.47(1) 




12.02(2) 




11.17(4) 


10.92(2) 


early-B 


5 


16:17:02.63 


-50:46:56.6 


14.43(3) 


14.26(8) 


12.78(3) 


12.51(7) 


11.64(4) 


11.17(7) 


early-B 


6 


16:17:03.7 


-50:47:05 


14.61(3) 
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Fig. 1. — Magnitudes comparative diagram Mjhk{2MASS) x MjHK{CamIV). The continuous line shows 
the relation expected if the two photometric systems were equal. 
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Fig. 2. — Combined false-image made from the J. H and rihK LNA images. North is to the top east to the 
right. At the botton are detailed views of the infrared nebulae J [le ft) [center) and nhK{right) images. 
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Fig. 3. — nbK band contour map of the infrared nebulae associated with IRAS 16132-5039. The contours 
start at 2.2x10"^ Jy/beam, with the same intervals (the beam size is 2 x 2 pixels). The positions of selected 
infrared sources refered in Table 2 and the IRAS coordinate elipse error (dotted line) also are indicated. 
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Fig. 4. — Contour diagram from the A MSX band image (8.28/Ltm), overlaying a LNA's H band image. The 
contours start at 2.8x10"^ W/m^sr and are spaced by 8xl0~^ W/m^sr. 
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Fig. 5. — Spectral energy distribution of the "A" and "B" MSX sources. The mid infrared data (squares) 
were obtained from the integrated fltix of the individual sources in the MSX images (bands A=8.28/xm, 
C=12.13;um, D=14.65;um and E=21.34/im) while the far infrared data (triangles) were taken from IRAS, 
and are listed in Table 1. 
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Fig. 6. Diagram of the spacial distribution of sources detected at LNA's H band image. The two regions 
(a) (cluster) and (b) (control) are delimitated by boxes . 
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H-K H-K 

Fig. 7. — Color-color diagrams for two regions in our survey. The region labeled "control" (left panel) 
contains only foreground objects; the other labeled "cluster" (right panel) has also a foreground population 
(open circles) but presents objects with infrared "excess" (filled triangles). All objects with "excess" are 
located in the " cluster area" . The locus of the main sequence and giants branch arc shown by the continuous 
lines taken from Koornnecf (1983). while the dashed and dotted lines follow the reddening vectors taken 
from Rieke & Lebofsky (1985). The location (plus signs) of Ay — 0, 10, 20 magnitudes of visual extinction 
also are indicated. The cluster members candidates are labeled by numbers. 
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Fig. 8. — The J versus (J — H) color-magnitude diagram of the sources in table 2. The locus of the main 
sequence at 3.7 kpc is shown by the continuous line. The intrinsic colors were taken from Koornneef (1983) 
while the absolute J magnitudes were calculated from the absolute visual luminosity for ZAMS taken from 
Hanson et al. (1997). The reddening vector for a BO ZAMS star (dotted line) was taken from Rieke & 
Lebofsky (1985). We also indicated the location (bold numbers) of Ay = 10 and 20 magnitudes of visual 
extinction as well as the sources that show "excess" (open up triangles) and do not (filled down triangles) 
in the color- color diagram. 
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Fig. 9. — The J, H and nhK contour maps of the infrare nebula region. The contours start at 0.44 (J), 
{H) and 2.2xl0~* Jy/beam (K), and are spaced by 0.37, 1.8 and 1.1x10"^ Jy/beam respectively. 
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